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The invention relates to a method and a unit for substracting quantization 

noise from a PCM signal. 

Pulse code modulation (PCM) is a widely used technique to represent e.g. 
audio signals in the digital domain. Linear PCM uses evenly distributed quantization levels 
5 whereas the logarithmic PCM uses a distribution where small-signal amplitudes are 
quantized using more closely spaced levels than large-signal amplitudes. An example to 
linear PCM is the CD-audio, which uses 16 bits to represent uniformly distributed signal 
levels. 

Fig. 5 illustrates the quantization of a linear coded PCM signal amplitudes 
10 (curve (a)) and of two logarithmically coded PCM signal amplitudes (curve (b) and (c)). 
Cuve b illustrates the quantization of ^l-law coded PCM signal amplitudes, wherein the M 
has an input-output magnitude characteristic of the form: 



log(l + (1) 



15 I U 

log(l + A^) 

where \s\ is the magnitude of the input, \y\ is the magnitude of the output and /i is a parameter 
that is selected to give the desired compression characteristic. 
20 In Fig. 5 curve c) illustrates the quantization of A-law coded PCM signal 

amplitudes, wherein the A-law has an input-output magnitude characteristic of the form: 



I J =111^, ^ (2) 

'•^ 1+log^ 

wherein the meaning of y and s is as mentioned above and A is constant. The prior art ^-law 
and A-law coders are standardized by the Comite Consultatif International de T616fonique et 
' Tel6graphique CCITT as the G.71 1 64kbit/s PCM. The quantizer characteristics of the G.71 1 
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are disclosed in-. N.S.Jaya„.andP.Noll.-T>igital coding ofWaveto^ns^E^glewoodCli^^ 

NJ- Prentice-Hall, 1984. . 

Thequantizedsamplevaluesofasignalmaybemodeledmathematxca^^^ 

S*[w] = 5W + g[w], 

5 

whereinsnnlreprcsents .he qua„««dvalueofanoriginalun,uan.ized sample val^^^ 
L<^„;iL.sthe,ua„.iza«o„e™r,esul«ng.o.said,— .The,^^^^^ 
eoor can be interpreted as undesired additive „oisede™fl.c quality offl.eonga«l 

10 sample values s[n]. 

Thus, it is the object of the invention to provide a method and a unit for 
substracting said quantization noise ftom a PCM signal. 

This object is solved by a method comprising the steps of: 
calculating for each fame of said PCM signalaquantizadonnoise level B, according to the 

15 following equation: 



20 



wherein 



S*min[n] 
25 S*max[n] 



• indicates a specific sample of the PCM signal; 

: represents the minimum quantization noise level for a specific sample value 
s*[n]ofsaid PCM signal; 
: represents the maximum quantization noise level for the specific sample 
value s*[n] of the PCM signal; 
w[n] •• represents a window-function; and 

W : represents the number of samples per window; 

30 and 

substracting the quantization noise as representedby said quantization noise level Bq ftom 
said PCM signal. 
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It is the essential idea of the invention to calculate said quantization noise 
levels Bq for each frame according to the claimed formula and to use any suitable 
(background) noise substracting method known in the art - perhaps after carrying out slight 
modifications - to substract the quantization noise from the PCM signal. 
5 Throughout the whole description * indicates the content of quantization noise 

in the co-ordinated quantity or signal. 

The window-ftmction may be chosen arbifrarily. 

The claimed formula requires the minimum quantization level S*min[n] as 
well as on the maximum quantization level S*n.x[n] of the sample values s*[n] of the PCM 
10 signal. 

According to a first advantageous embodiment of the invention the minimum 
quantization level S*„^„[n] as well as the maximum quantization level S*™«(n] are known 
from the quantization method applied for generating the PCM signal. 

According to a second embodiment said quantization levels STnunM and 
15 S*^[n] are not known and must therefore be predicted according to the formulas claimed in 
claim 3. 

Advantageously neither in the first nor in the second embodiment any 
additional information must transmitted together with the PCM sample values snn]toena^^^ 

the conduction of the claimed method. This is because in the first embodiment the 
20 quantization levels and in the second embodiment reference levels i can be assumed to be 

known (see the last two paragraphs above) when the method for substracting the quantization 
noise according to the invention is carried out. Because of that advantageously bandwidth is 
saved. 

Preferably, the noise substraction is based on a background noise substracting 
25 method and system as described in the non-prepublished European Patent application having 
the application number EP 01201304.1 of the same appUcant. The application of said system 
is described later as preferred embodiment. The whole content of said patent application shall 
be regarded as part of the disclosure of the present application. 

Finally, the above-identified object is solved by a noise substracting unit 
30 according to the subject matter of claim 9. The advantages of said unit correspond to the 
advantages mentioned above refening to the clauned method. 

Preferably, the claimed unit and method are applied at the decoder side, i.e. 
just before the PCM signal is decoded again. However, nevertheless it might also be appUed 
at the encoder side. 
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The description is accompanied by five figures, wherein 
Fig. 1 showsaquantizationnoisesubstractingunitaccordingtotheinvention 

including the background noise subtracting unit; 
Fig. 2 shows an example for illustrating the correlation between reference values I 

and sample values S*[n]; 
Fig. 3 shoes a detailed illustration of a background noise substracting unit; 

Fig. 4. 1 . shows the spectrum of a clean male voice; 
Fig 4 2 shows the n-law coded version of the spectrum shown in Fig. 4.1.; 
Fig. 4.3. shows the ^-law coded version of the spectrum shown in Fig. 4.1. after 

quantization noise substraction according to the invention; 
Fig. 4.4. shows the usage array u[i] generated for applying the quantization noise 
substraction method to the spectrum shown in Fig. 4.2; 
showing quantization levels for linearly and non-linearly coded PCM-signals 

as known in the art. 

In the following apreferred embodiment of the invention v«ll be described by 

referring to Figures 1-4. ^ 

Fig. 1 showsa,u»ntization noise substracting unitlOO-VOOforsubstraetmg 

20 quantizationnoiseftomaninputlinearlyorlogarithmicallycodedPCMsignal. 

First of all said unit comprises a quantization noise level calculatmg umt 100 
forcalculatingaquanti.ation„oiselevelBqforeachframeofsaidPCMsignalunder 
consideration ofbothaminimumandamaximum quantization noise level S«^„[n).S ™.lnj. 
More specificaUy.the quantization noise levelBq is constant for each tone of said PCM 

25 si^ and calculated individually for a specific frame according to ti>e following equation: 



Fig. 5 



wherein 
n 

S*min[n] 



• indicates a specific sample of the PCM signal; 

: represents the minimum quantization noise level for a specific sample value 
s*[n] of said PCM signal; 



S*^x[n] : represents the maximum quantization noise level for the specific sample 
value s*[n] of the PCM signal; 
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w[n] : represents a window-function; and 

W : represents the number of samples per window comprising at least the specific 

frame for which Bq is currently calculated. 

5 Typically.thenumberofsamplesperwindowistwicethenumberofsamplesper 

\js^n]-S^Jn])'W[n]Y_^ represents the mean 
In equation 4 the term ^ 2 

.quare value ofthe quantization error. i.e. the energyoffte quantisation noisein .he cunen^ 
frame. 

The frame size is typically 2ms. 

According to a first embodiment said quantization noise levels S*„.„[n]. • 

S*max[nl are assumed to be known. 

However, according to a second embodiment of the invention these 
quantizationlevelsSWn].S*^[n]are not assumed to be la.own.e.g. because theori^^^^ 
15 signalsmighthavebeenquantizedwithdifferentquantizersandtheremayhavebeenextra 
scalings applied on the samples before or after coding. In these cases the quantization 1^^^^^ 

have to be predicted. 

To predict the quantization levels, tire entire original signal, heremafter also 
„ferredtoaswavet6m,isscanncdtoreceivethesamplevaluessnn] of the quantized 
20 wavefom,.Subsequenfly,acheokismade.ifallpossiblerepresentationlevelsipredefined 
due ti,e quantization method appUed to the waveform arerepresentedCused) by at le^to^^^^ 
saidsamplevaluessnnlFora 16-bitlinearPCMquantizationmett«,dfl.ereare2 =65536 
possible representation levelsiandflie result of said cheek is described inaso called usage 
tableu[i]with65536values.Thevaluesofsaidusagetableu[i]areeitherl(itthere.sa. 
25 leastonesamplevalueofthewholewaveformcorrespondingtothecurrentrepresentation 
level i) or 0 (if there is no sample value of the whole waveform corresponding to the current 
representation level i). Expressed mathematically 



30 r-1 „,Ji^P' -r-'<i<r-' (5) 

«['l = ™\'-ll, otherwise) 
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In formula (5) i is the wavefonn length, i.e. the number of samples of the 
whole waveform, and N is the number of bits used by the Unear PCM quantization; in the 
above example N= 1 6. This operation can be simplified by initializing all values oIu[i] wrth 
0 and using fonnula (6): 



0 < n < L 



(6) 



10 Now equations (5) or (6) are used to find special representation levels i^M 

^^^i^Mrequiredforpredictinsthedesiredquant^^^^^^^ 

isdefmedasthatrepresentationlevelwhichis-startetfromiW-thenextsmallernon-z^ 
representation level for which u[n]=l and i^[n]is defined as that representation levd 
is - started from i[n] - the next bigger non-zero representation level for which u[n]-l . 
^5 With the help of said co-ordinated representation levels the desired 

quantization levels are predicted or calculated according to the following equations: 



s*min = i[n]-(i[nl-imin[n])/2 
S'™uc = i[nl + (imax[n]-i[n])/2 



(7) 
(8) 



25 



wherein . , , 

i : represents one out of a plurality of predefined possible representation levels; 

i[n] : represents that predefined representation level which corresponds to the 

sample value s*[n]. 
Fig. 2 illustrates the calculation of the usage array, of said special 
representation levels i^M, and of the desired quantization levels S*„.„[n]. S Wn] 

for a given wavefonn consisting of L=9 sample values S*[n] which are in the example 
distributed as shown in Tab. 1 : 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


s*[n] 


3 


2 


-1 


-1 


2 


1 


I 


3 


-3 


=i[n] 





















Tab. 1 
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In Fig. 2 there are provided 8 different possible representation levels i= -4, 
1, 0, 1,..., 3 resulting from a N=3 bit linear PCM quantization. The usage array u[i] calculated 
according to equation (5) or (6) for the example of Tab. 1 is: 

u[-3] = u[-l] = u[l] = u[2] = u[3] = 1; and 
u[.4] = u[-2] = u[G] = 0 

Now, from said usage array the reference levels iminW and imaxW are 
calculated e.g. for n=2 with s*[2H[2] =1 - according to the definition given above - to: 

imax[2]=2andimin[2] = -l 

and the desired quantization levels result -according to equations (7) and (8) to: 

s*.d„=l-(l-(-l)V2 = 0 
S'max=l+(2 -1)/2 = 1.5 

Special care has to be taken when i[n] co-ordinated to s*[n] is the first or last 
non-zero value in the usage-array u. One option is to allocate two extra array values uf-f'''- 
l] and set them both to 1 and always continue the search until a non-zero value is 

found. 

Another question is what to do when some valid quantization levels aren't 
used in the waveform. Especially waveforms shorter than 1 second will result in a usage- 
array with many unused (but valid) quantization levels. Since this usually will give a wrong 
result for the quantization step-size at boundaries (when the waveform is at its maximum or 
minimum value), the effect on the measured level of quantization noise will be small. 
Although an algorithm could be used to predict the unused valid quantization levels, the 
added value of the extra complexity is not worth the performance gain, m practice, it seems 
to be enough to assume that the used signal levels in the waveform are the only valid levels. 

After the quantization noise level Bq has been calculated either with the 
known or with the predicted quantization levels s*™„ and s*„^ it is input into a background 
noise subsfracting unit BNS 200 which is also part of said quantization noise substracting 
unit. 
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Fig. 3 Shows a detailed illustration of said background noise substracting unit 
200 which is embodied to substract the quantization noise from a PCM signal in the 
frequency domain. It comprises a magnitude forming unit 210 for forming the magnitude 
|S*[k]| of the spectrum sm of the PCM signal, a signal-to-noise ratio calculating unit 220 
for calculating the signal-to-noise ratio of the PCM signal according to the formula: 



SNR[k] = |S*MPq, 



(9) 



wherein Bq is noise level representing the level of quantization noise in said PCM signal. 

Said signal-to-noise ratio SNR[k] is input into a filter update unit 230 for 
calculating a filter magnitude F[k] according to a predefined filter algorithm based on at least 
one filter update parameter. Said filter magnitude F[k] is input into a calculation unit 240 for 
calculating an output spectrum by multiplying both the real part R{S*[k]} and the 
imaginary part I{S*[k]} of the spectrum S*[k] with said filter magnitude F{k}. The filter 
update parameter and thus also the filter magnitude F[k] are preferably adjusted such that 
output spectrum S^M is at least substantially free of quantization noise. 

Because said background noise substracting unit 200 is embodied to carry out 
the substraction of the quantization noise in the frequency domain, the quantization noise 
substracting unit fiirther comprises several units for computing the spectrum S^k] of the 
PCM signal to be input into said backgromid noise substracting unit 200 and for transforming 
the spectrum S^[k] output by said unit 200 back into a time domain signal. These units are 
described in the following paragraphs. 

Referring back to Fig. 1 the quantization noise substracting unit fiirther 
comprises a first window block 300 for generating a weighted PCM signal s*w[n] by 
weighting the frames of the input PCM sgnal with a fu^t window wl[n]. Said weighted PCM 
signal s*w[n] is input into a spectrum computing unit 400 for computing the spectrum S*[k] 
from said weighted signal. Preferably, the spectrum is computed by using the Fast Fourier 
Transformation algorithm. The spectrum S*[k] is than input into the background noise 
substracting unit 200 and processed therein as described above by referring to Fig. 3. 

The spectrum S''[k] output by said background substracting unit 200 is input 
into a re-transfoiming unit 500 for being transformed back into a signal s^M in the time 
domain. In the case that the spectrum S*[k] has been calculated by using FFT the re- 
transformation is done by using the inverse FFT. 
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The quantization noise substracting unit may further comprise a second 
window block 600 for generating a weighted output signal s^[n] by weighting the signal 
s^[n]received after the transfomiationwithasecond window W2[n] and an overlap and add 

unit 700 for individually calculating a final output signal ^ [w] for each frame of the 
PCM-signal from said weighted output signal sVn] such that the transition between two 

successive output frames is smoothed. The first and the second window Wl and W2 are 
preferably identical. They are independentofthe window-function w[n] used for calculatmg 
the quantization noise level Bq.However.it is preferred to use the same fimction for w[n] as 

forWl. 

The provision of the overlap and add unit 700 only improves the quality of the 
output signal but it is not absolutely necessary for carrying out the quantization noise 
substraction;consequently.itmight be leftoutTHe same applies to the first and the second 

window block 300 and 600 which are both only provided to improve the operation said 

overlap and add unit 700. 

Fig. 4 shows an example illustrating the good signal quality achieved when 

applying the claimed method to a ^i-law coded speech signal. 

More specifically, Fig. 4.1 shows the spectrum of a clean, i.e. free of 
quantization noise, malevoice recording wUh the utterance:"Any dictionary will giveat least 

twenty different definitions for it". Fig. 4.2. shows the ^-law coded version of the spectrum 
shown in Fig. 4.1. From that Fig. 4.2. it can be seen that said ^-law coding causes 
quantization noise added to theclean spectrum ofFig.4.1.In Fig. 4.2. the quantization noise 

is represented by the grey-shaded areas 410; the corresponding areas in Fig. 4.1 are clean. 
Fig 4 3 shows the spectmm from Fig. 4.2 after the quantization noise substracting method 
according to the invention has been applied. It can been seen that this spectrum substantially 

corresponds to original clean spectrum shown in Fig. 4.1 ; that means that the method 
according to the invention is well able to clean the spectrum of a PCM signal from 
quantization noise. Fig. 4.4shows the usage array u[i] as generated bythe method according 
to the present invention for reducing quantization noise from the spectrum shownmFig. 4.2. 



